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Starting from the Liouville-von Neumann equation, under a weak coupling limit we derive the 
Lindblad master equation for the one-dimensional quantum Ising model in a Markov approximation 
and a rotating wave approximation. We also prove that the steady solution of the Lindblad equation 
is the canonical distribution independent of the dissipation rate. 

I. INTRODUCTION 

The Lindblad master equation plays an important role in open quantum systems. As an example, it has been 
widely used in quantum optics^. It was originally derived by Lindblad using quantum dynamical semigroups^. Attal 
and Joye^ obtained it through taking a continuous limit of repeated interactions with a sequence of baths in a 
given density matrix state. Recently, Brasil, Fanchini, and Napolitano^ provided a relatively simple derivation of 
the Lindblad equation starting from a general Hamiltonian and the Liouville-von Neumann equation in a Markov 
approximation and a rotating wave approximation under weak-coupling limit. 

For many-body systems, many studies focus on the Lindblad equation for open systems interacting with the envi- 
ronment at their boundaries -^. Here we derive the Lindblad master equation for a one-dimensional quantum Ising 
system interacting with a heat bath at each site following the schema given in Ref. 1. 

Our derivation and discussion of the Lindblad master equation will proceed as follows: In Sec. [H] we diagonalize the 
Hamiltonian for a composite system consisting of a quantum Ising chain and a heat bath and write it in the interaction 
picture. In Sec. IIIII we start from the Liouville-von Neumann equation describing the evolution of the density matrix 
of the composite system to derive the Born-Markov equation under a weak system-environment interaction limit and 
a Markov approximation. These results are then combined to obtain the Lindblad equation for the quantum Ising 
model with the bosonic heat bath interacting at each site in Sec. IIVI The steady solution of Lindblad equation is 
shown to be the canonical distribution independent of the dissipation rate in Sec. [Vj followed by a brief summary in 
Sec. EH 

II. HAMILTONIAN AND ITS DIAGONALIZATION 

The Hamiltonian for the quantum Ising chain interacting with a bosonic heat bath is 

H = H S + H B +H SB , (1) 

with 

N N-l 

H s = -h x '£<T?-j'£a*a? +1 , (2) 

i=l i=l 

H B = ^u) l) tf l}i bf H , (3) 



H SB = Y / MbU + b^K, (4) 
p.i 

where H$, H B , and H$ B are the Hamiltonians for the quantum Ising chain, the heat bath, and their interaction, 
respectively, TV is the total number of spins, o\ represents the Pauli matrix along j direction at site i, h x is a traverse 
field, A/3 is the coupling strength, (6^) creates (annihilates) in mode j3 with an energy ujp a boson coupling to the 
spin at site i. Equations ([T]) to (|4|) have been studied in Ref. [7] using Green's functions. 
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In order to diagonalize Hs, we first apply a Jordan- Wigner transformation 8 to Hs and Hsb- After the transfor- 
mation, we have 

H s = -h x Ef =1 (l - CjCj) - JEf =1 (C]c] +1 + C]C 3+l - C,C] +1 - CjC j+1 ), (5) 

Hsb = ]T Wpi + W - C j c i)> ( 6 ) 
where Cj and Cj are creation and annihilation operators at site j for the Jordan- Wigner fermions, respectively. 



Next, a Fourier transformation leads to 



q>0 



-i smq — j- + cos g 

H ™ = S K^+o - ^+^1) + w<3U^ - . ( 8 ) 

ffn = 5>„&£ g &/» 9) (9) 

where C<j and C q are the corresponding operators in momentum space. 
Further, a Bogoliubov transformation^ of 



results in 



with 



and 



HjSB " VN h 



Vq = u* q C q + v_ q Cl q (10) 

9 

3 

H SB = ^2 H jSB , (12) 

= x p( b h a rti + 6 « a ]fe<?)' ( 13 ) 



aifc 9 = (u* k u k+q - v k+q vl)(r)lr}k+ q - Vk+qVk), 

aikq = (V-kUk+q - Vk+qU-kjV-kVk+q, 

a 3 k q = {u* k v*_ k _ q - u*_ k _ q vl)rilr)l k _ q , (14) 



where f2 9 = 2JyJ (h x /J) 2 + 1 — 2{h x / J)cos q. 
In the interaction picture defined by 

= e^+^Ae-^+v^, (15) 

for an operator A, 

Hsb (0 = 4EE A ^ ^-Oifcge-^' + &«at^ e *w) , (16) 



V i=i fc,g,/9 



where the superscript I denotes the interaction picture and uijkqfi = tUj kq — up with uj\ kq = fl k + q — Clk, ^>2kq = 
Q-k + &k+q, and uJskq — — Qk — Q-k-q- In Eq. (|16[) . we have made use of the relations 

exp(iHst)a jkq exp(-iHst) = a jkq exp(-iu> jkq t), exp(iHst)a' jkq exp(~iH s t) = a] fcg cxp(iu> jkq t), (17) 

derived from Eqs. (fTTj) and (IT4"1) . 
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III. DERIVATION OF BORN-MARKOV EQUATION 

The composite system evolves according to the Liouville-von Neumann equation 

^ = -i[H,p], (18) 

where p is the density matrix of the composite system. In the interaction picture, the Liouville-von Neumann equation 
becomes 

^- = -i[H I SB {t), /(*)]. (19) 
Integrating Eq. from to t and substituting the result back into it gives 



dp T (t) 
dt 



(20) 



where p 1 (0) = p£(0) ® p B {$). 

Our goal is the evolution of Pg(t) — Tr B [p (t)], the density matrix operator of the system itself. So, 



at 

Assuming 



-iTr B [F| s (t)y(0)] - T± B [H x SB (t), J dt'iH'ssit'^p^t')}]. (21) 



Pb(0) 



t r pxd ( " w ^ h ^ 

ll/3q eX Pl k B T 



Tr s 



EU exp 



-iupb^bp,, 



Pq^ y \ k B T 



(22) 



for an equilibrium bath at the temperature T (k B is Boltzmann's constant), one sees that — iTT B [Hg B (t), p 1 (0)] = 
because T^B(bp q p B (0)) = Tr B (bpqP B (0)) = 0. Also, in Born approximation of weak system-environment interaction 

limits, we can write p 1 (t) = p$ (t) (g) p B (0) as the reduced density matrix for the bath changes little with time because 
of its huge size. So, Eq. (|2ip becomes 



t 

2EM = - T r B [H I SB (t),J dt'lHi^lpUt')®^]], 



(23) 



t' t" 

Ps{t')-Ps{t)=- J dt'TrslHUt"), J dt'^H^it" 1 ),^" 1 )®^]], (24) 
t o 

after integration. We see that the difference between Pg(t') and Pg(t) is of second order in Hg B . As a result, in the 
weak system-environment interaction limit, we can replace p 1 (t') in Eq. (1231) with p 1 (t) and arrive at 

t 

^fi = -Tr^HUt), J dt'm B (t'),pUt)®p B ]}. (25) 
o 

In the Markov approximation, the system is memoryless. This can be obtained by changing the integrated variable 
to t — t' and sending the upper limit of the integral to The result is a Born-Markov equation, 

°ekp.=-Tr B [H' SB (t), [ dt'iHUt-t'), Ps(t)®P B }]- (26) 
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IV. DERIVATION OF LINDBLAD EQUATION 



In this section, we shall derive the Lindblad equation from the Born-Markov equation (|26[) in the rotating wave 
approximation^. 

We first consider the contribution from Hig B . Denoting the corresponding density matrix as substituting 
Eq. (HU) into Eq. ([26]). and using! 

Tr B (b^b 0w p T B ) = Tr B (bl q bl >q ,p B ) = 0, (27) 
with (71(10/3)) = (bp q bp q ) = l/[exp(u;p/k B T) - 1], we have 

— .X. 

N 

x { [aL 9 aifcwP(s((«(^)> + l^wt-iuxvvlt-e) + a lfeg a| fc , 9 p{ s (n(^))e- <w «^*+^^'^(*- t ') 

« t lfcg p(s«ifc' 9 («(^))e ,; — ^-"i^C-O + a lfeg p{ s al fc , g ((n(^)) + i) e — ^+—^ C*-0 
- at fc ,^ s a lfc9 (n(^))e--«^ t +--'^( t -*') + au^a^n^)) + i) e ™<— '„C*-0 
+ [pi^L^^C^M) + l)e— ^ t+ -«'^( 4 -*') + /9 ( saifc , g a t lfc? (n( W/3 ))e^^ t — '«*('-*')] } . (28 ) 

To proceed, note first that the integral 

00 

/ dt'e^x'^' = nS(uj lklq p) ± iV—-— (29) 



dpi, 
dt 



-\-00 



kk'q/3 



indicates that only the values of fe', 9, and j3 satisfying uiik'qf) = contribute to the real part of the righthand side of 
(|28p due to the (5 function. Further, after the integration, one is left with exponentials of e l ( u ' lk 'if 1 ~ u ' lk '> f> ' t , which, in 
the rotating wave approximation 10 , vanishes unless uik'qp = uikq/3 or fc = k' . So, Eq. (j28|l becomes 



= - i[H 1LS ,p{ s } - [likq((n(uJik q )) + l)(P5 a Ife 9 «ife<? + a lkq a ikqPs ~ 2ai/c 9 Ps a ifc^ 

kq 



(30) 



where H 1LS = (1/iV) EfcgC^ifeg^ifeg' a ife<?] + Aw i/cg a u g a ifc<?) with ^1*9 = XL T >c u( n {^)) l( u ikq ~ Aw lfc , = 
^Pcu/fwii:, - w), and c w = X^ A|£ W , W(3 , 7ifc<7 = ^n,/^). and c = J2u, c <^/ N ' ( N ' is tne tota l number of 
boson states), which is a dissipation rate describing how fast the system dissipates to equilibrium. 

Similarly, the term containing only -ff^ss contributes an equation similar to (|30p because of the modes satisfying 
U2kq/3 = W2fc'ij/3 and again k = k'. The cross terms from H[ SB and ff^ss should vanish in the same approximation 
because the combination of fc', k and q to meet uJikqp = ^>2k'qp or u!2kqp = ^lk'qp, if any, is much less than that 
satisfying either utikqp = wifc' 9( 3 or Ld2kqp = ^2k'qp- Terms containing H^ SB do not contribute to the real part of the 
righthand side of Eq. (|28|) because uj^kqp is always negative. 

Collecting the relevant terms, we then obtain the Lindblad master equation in the interaction picture 



M 

dt 



A 

C Y1 Yl ["fjkq((n(Ujkq)) + l){Psa]kq a ikq + Ojkg^kqPS ~ 2a jkqPW 3 kq) 



3 = 1 k,q 
2 



(31) 



: E5Z \jjkq((n(u jkq ))(p I s a : j kq a] kq + a^a]^ - 2a) kq p s a jhq ) 



3 = 1 k,q 
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where k and q take values satisfying ujjkqp — 0. In Eq. (|31j) . we have neglected a Lamb shift term —i[Hr,Sip] with 
= H\ls + ^2LS + ^3LS because it is of higher order—. Transforming back to the Schrodinger picture, we find 

dp 2 

-JjT =- *[^s,Ps] - C EE \ljk q {{n{ujk q )) + l)(p s a] kq ajk q + a] kq a jkqPs - 2a. ]kqPs a] kq ) 

(32) 



2 

C 

3=1 fe,<3 



[^?(M w ifc?))(Psajfc g aj- fc g + a jkq a] kq p s ~ 2a] kqPs a ]kq ) 



In order to write Eq. (|32p in the familiar Lindblad form, for two energy levels E\ and E m with £"„> — E\ = Wjfc g , 
we may let Vm^z = Ojfeg and V^_^ = which are respectively thermal jump operators representing emitting and 
absorbing a particle with energy ujj kq and jumping to a lower and higher energy state. Also, let Wi^ m = 7jfc 9 (^(wjfc g )) 
and W m -^i = Jjkq((n(ujj kq )) + 1), which are transition probabilities from Zth state to the mth and vice versa, so thatjii 

W(»K fcg ))) _ exp (za^U exp f^LZ^V (33) 



W m _>.j 7jkq({n{u)jkq)) + 1) \ k B T J \ k B T 

Equation (j3"2")l then becomes 

^ = -i[ff s , P5 ] - c E W^ m {V^ m V^ m p s + Ps V^ m V^ m 2V^ m p s V l lj 

rn > I 

-¥rn) (34) 

= -i[H s , ps]-c W^ m (V^ m V^ m ps + PsV^ m V^ m - 2V l ^ mPs V l \ m ), 

the usual Lindblad form. Although the form of the transition probability Wi^„, depends on the environment and 
determines the details of the process, universal properties only rely on Wi_> m /W m ^ir^- For example, in Ref. (12), 
Wi^ m — /3 m with /3 m the probability for the system to stay in the mth state in equilibrium. This completes our 
derivation of the Lindblad equation for the quantum Ising model. 



V. CANONICAL DISTRIBUTION IS THE STEADY SOLUTION OF THE LINDBLAD EQUATION 

In this section, we discuss the steady solution of Lindblad equation. 

The meaning of the Lindblad equation (|34l) is clear. If the second term in the right hand side is neglected, it is 
the quantum Liouville equation determining the quantum fluctuations of the evolution of density operator p$; while 
if the first term on the right hand side is neglected, the diagonal part is 

^ = c^Wj H p j: j ~ Wi^jpu), (35) 
which is the classical master equation, and the off-diagonal part is 

% 1 = -|(E^ + E^)P«> (36) 

kjti Ift 



which decays exponentially. Thus the Lindblad equation (|34|) naturally integrates the quantum and thermal fluctua- 
tions together. 

It can be readily checked that for a time-independent Hamiltonian, the equilibrium density matrix of the canonical 
distribution pe = exp(— iJs/fcsT)/Tr[exp(— Hs/k B T)] is the steady solution of the Lindblad equation^. To this end, 
we note first that there is a detailed balance condition, 

m v it u exp(-H s /k B T) ( H s \ ( H s \ 

W^ mPE V m ^ = ljk qH^kq)) ^ [eM _ Hs/kBT)] W«P [j-f ) «P {-J-f ) 

, t « f"jkq\ exp{-H s /k B T) (37) 

= j jkq {n{u) jk q))eicp ) a i k i' 



'Tr[exp(-H s /k B T)] 



Wm^lVm^lPE, 
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where uses have been made of Eqs. (1TH) and (j33l) . Then, on the righthand side of Eq. (|34|). the first term [Hs, pg] = 0. 
For the second term, using the property of the thermal jump matrices V;_> m = V^L^ from their definitions and 
substituting Eq. (1571) into the righthand side of the Lindblad equation, the second term can be explicitly written as^, 

Wl^ m {V^ m Vl^ mP E + PEV^ m Vl^ m - m^mPEVUra) = Wl^ m V^ m Vl^ m pB ~ Wl-> m Vl^ m V^ m p E , (38) 

and 

Wm^liVl^Vm^lpE + PEVl^V m ^i - 2V m ^l PE V^ l ) = W m ^iVl^V m ^ipE ~ W m^lVm^iV^p E ■ (39) 

These two terms cancel with each other and thus the second term of the right hand side of the Lindblad equation (|34|) 
equals zero too. This derivation is similar to that from Ref. (3), but we take the thermal jump involving all energy 
levels into account. Thus for the weak coupling situation, the Lindblad equation reduces to the canonical distribution 
solution in long time. Note that this steady solution does not depends on c. 

VI. SUMMARY 

We have derived the Lindblad equation for the quantum Ising chain weakly interacting with a heat bath. Further 
we have confirmed that the steady solution of this equation is the equilibrium canonical distribution independent of 
the dissipation rate. 
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